On the Eigenvalue Density of the non-Hermitian Wilson Dirac Operator 
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We find the lattice spacing dependence of the eigenvalue density of the non-Hermitian Wilson 
Dirac operator in the e-domain. The starting point is the joint probability density of the corre- 
sponding random matrix theory. In addition to the density of the complex eigenvalues we also 
obtain the density of the real eigenvalues separately for positive and negative chiralities as well as 
an explicit analytical expression for the number of additional real modes. 
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Introduction. In the past two decades, there has been 
an increasing interest in non-Hermitian random matrix 
theory (RMT) [1]. To name a few apphcations, quantum 
chaos in open systems [21, dissipative systems [3| and 
QCD at finite chemical potential Some features of 
the model we are considering also occur in the condensed 
matter system analyzed in Ref. Q. 

The connection between the infrared limit of QCD and 
RMT has been well understood in the continuum limit 
since the early 90 's [6]. It is based on the universality 
of chiral RMT in the microscopic limit (or e domain) [7| 
with chiral RMT described by the same chiral Lagrangian 
as QCD. The main advantage of RMT is the availability 
of powerful methods to derive analytical results, and re- 
cently this approach was applied to QCD at finite lattice 
spacing |8l4ll] . It was shown that the e limit of the chiral 



Lagrangian for the Wilson Dirac operator Dw tl3|, |14| 
can be obtained from an equivalent RMT. Discretization 
effects of the spectrum of Dw have been studied directly 
by means of chiral Lagrangians [lij . but using 

RMT methods will enable us to obtain results that were 
not accessible previously. 

The aim of this paper is to obtain analytical results 
for the eigenvalue density of Dw for the RMT model 
proposed in Ref. We consider the quenched case. 

RMT. We consider the random matrix theory 



aA W 
-W'^ aB 



(1) 



distributed by 



P{Dw) oc exp — — (trA^ + trS^) - ntiWW'^ . (2) 

The matrices A and B are Hermitian n x n and {n + 
u) y. in + v) matrices, respectively, and the entries of W 
are complex and independent. In the microscopic limit, 
with rt — >■ (X) at fixed rescaled eigenvalues z = 2nz and 
lattice spacing = na'^/2, the spectral properties of 
this RMT become universal and agree with Wilson chiral 
perturbation theory in the same limit (with n identified 
as the volume of space-time) apart from the squared trace 
terms 14, 1^. The finite integer < n is the index of 
the Dirac operator and is kept fixed. 



The matrix Z?w is 75 = diag(l„, — l„+,y)-Hermitian, 
i.e. = 75-Dw75- Therefore its eigenvalues are either 
real or come in complex conjugate pairs. The v generic 
zero modes at a = become the generic real modes of 
-Dw at finite lattice spacing. Furthermore, Dyj may have 
2n — 21 additional real eigenvalues which appear when a 
pair of complex eigenvalues collides with the real axis. 

In Refs. [9l4l2| the technically simpler case of the Her- 
mitian Wilson Dirac operator D5 = 75DW was studied. 
Although sgectra of D5 have been studied in the lattice 
literature [17|, only the eigenvalues of Dw are directly 
related to chiral symmetry breaking which is our main 
motivation to study its spectral properties. 

The joint probability distribution (jpd). To preserve 
the 75-Hermiticity of Dw we can only quasi-diagonalize 
by a non-compact unitary matrix U G U(n, n + v), 
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In contrast to the diagonalization of a Hermitian matrix 
such as D5, the matrix X may only be quasi-diagonal 
where xi, X2, 1/2 and X3 are diagonal matrices of di- 
mension n — I, I, I and n — I + v with < I < n the 
number of complex conjugate pairs. The complex eigen- 
values are given by iz,z*) — {x2 + iy2,X2 — ^2)- The 
ensemble Dw decomposes into n-f-l disjoint sets of quasi- 
diagonal matrices ([3]) with a fixed number of real eigen- 
values. The joint probability density of the 2n-\-v eigen- 
values Z — (zii-, . . . , Z„r, Zi\, . . . , Zn+jyj) G C2n+iy y^Q 

obtained by integrating over U . This calculation will be 
discussed in detail elsewhere. We only give the result for 
1/ > which is not a restriction because of the symmetry 
V — )■ —V. The jpd is given by 



p{Z) CX A2n+v{Z) dct 
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FIG. 1: The projection of pc onto the imaginary axis for f — 
1 and a = l/\/8- The Monte Carlo simulation (histogram, 
bin size=0.4) contains 200000 matrices with n = 50. This 
simulation nicely confirms our analytical result (blue curve) 
and shows the deviations from the a = result (red curve) . 

where erfc is the complementary error function and 
6^^\x + iy) = 6{x)S{y). Due to 75 the permutation group 
&{2n + I') is broken to &{n) x &{n + ly) which reflects 
itself in the product of the Vandermonde determinant 
A2n+u{Z) and the other determinant in Eq. (U). The ex- 
pansion of the delta functions yields the jpd for each of 
the n + 1 subsets with a fixed number of complex eigen- 
value pairs. The two-point distribution 172 splits into one 
term for the real eigenvalues (72r and one for the com- 
plex conjugated pairs g2c as it is also known for the real 
Ginibre ensemble and its chiral counterpart JJ,]. 



The eigenvalue densities for the real and complex 
eigenvalues can be obtained by integrating over all eigen- 
values except one. The spectral density can be de- 
composed into the density of real modes, for posi- 
tive chirality ((V'ItsIV') > 0)j Pi for negative chirality 
((V'ItsI''/') < 0), and the density of complex pairs, pc, 



j P{Z) n ^[^j] = P\{xi\)5{yi\) 



Pc(zir) 



Pc{zi\) 



, (5) 



(6) 



Note that the chirality reflects the conventions of the 
RMT. By expanding the first row of the determinant 
in Eq. ^ and re-expressing the additional factors from 
A2n+iy(^) as A'f — 2 partition functions we obtain 

Pc{z) = g2c(z){z - z*)Z'^^^2{z,z*;a), (7) 
pA^) = J g2r{x,x'){x ~ x')Z'^^^2ix^x''^'^)dx' ■ (8) 

A similar factorized structure was found in Ref. [l9j . 

Expanding the first column of the determinant and in- 
tegrating over all eigenvalues except zu , we find the same 
expression for pc and the density pi of the real modes 
originating from g2 (using Eq. (6)). However, there is an 
additional contribution to pr due to the last v rows which 
gives the distribution of chirality over the real eigenvalues 



Px^ Pi- Pr- 



(9) 



Additional rows of some of the determinants have to be 
expanded to express them into known partition functions. 
We have checked for i/ — 1 and i' = 2 that the result 
agrees with previously derived expressions 



Hi. 



In the microscopic limit, the partition functions in 
Eqs. ([7]) and ([5]) can be expressed in terms of integrals 
over U(2). They can be simplified using the eigenvalues 
of the U(2)-matrices as integration variables resulting in 



= "",0 ^5/20- / e-(^°^'^^+^°^'^^)-''''(^°^''^^+™^''^^^sinc[y(cos(^i-cosvP2)]cos^.(^i+^2)i?^fc, (10) 

\2n/ (27r)^/^2a J[o,27r]2 

fx\ If exp[A?-A2]erf [Ai,V2Ai] -exp[A2-A2]erf [A2,v^A2] 

pA^I = -cosiy{(pi+(p2)D(fk, 11 

\2n/ IbTT^ J[Q 27r]2 COS — COS (^2 

Px(7r) = TTT-^ 1 s'{[siK^+i{si)I^{iS2) + iS2K^{si)I^+i{iS2)]— ——dsids2{l2) 

\2n/ loTTo^ Jr2 Si — ZS2 (^^ ^ 1)' 
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The functions sine, erf, //, Ki and J^'-* are the sinus car- dinalis, the generalized incomplete error function 
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FIG. 2: Along the imaginary axis the difTerence of pc for 
different v is much clearer at small a (solid curves) than at 
large a (dashed curves) where they are almost the same. 



(erf (6, c) — erf(c) — erf (6)), modified Bessel function of 
the first and second kind and the /-th derivative of the 
Dirac delta function, respectively. The integration mea- 
sure is induced by the invariant U(2) measure, DLpk — 
sin^(((pi — ip2)/2)d(pid(p2 and — 2a(^cos(pj — x/Sa^). 
Because of the (5-function only the algebraic singular part 
of the contributes to (which was already obtained 
in Refs. 0, [ll|). The distribution vanishes for j/ = 
and can be obtained from the generating function for the 
eigenvalue density oi ^^{Dw + 1^). [9| Comparisons of the 
analytical results with simulations of the random matrix 
model ([!]) are shown in Figs. [T] and ID The normalizations 
are chosen such that the integral over is equal to v. 
The other constants are already fixed by this choice. 

For small increasing a the complex eigenvalues move 
parallel to the real axis according to a Gaussian distribu- 
tion with a width of 2a (See Fig. [3]). Therefore the density 
of the projection of these eigenvalues on the imaginary 
axis is very close to the a — result. For large a the dis- 
tribution of the real parts of the complex eigenvalues de- 
velops a box- like shape from —8cP to 8a^ which can be de- 
rived from a saddle point approximation of Eq. (|10p (See 
Fig.[3|) and the oscillations disappear (See Fig. [2]). Along 
the imaginary axis pc{a ^ 1) becomes a~^erf (y/4a). 

Near the real axis pc behaves as y'^'^^ for small a but 
is linear in y for a large enough (See Fig. (2). In the 
continuum limit it peaks around the imaginary axis and 
eventually gets the form of the continuum microscopic 
eigenvalue density. 

Real modes. For a Wilson Dirac operator with index 
h' there are at least v real modes. The additional real 
modes result when complex conjugate eigenvalue pairs 
enter the real axis. The average number of these modes 
follows from the integral 
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FIG. 3: The distribution pc along a parallel axis to the 
X-axis (here aX y — 4Qa?) is Gaussian shaped for small o 
(solid curves) and develops a plateau for large lattice spacing 
(dashed curves). 



In the limits for small and large lattice spacing we find 



J2(iy+l)^ 

a, 



a< 1, 
a> 1. 



(14) 



This is shown in Fig. 21 For large lattice spacing the 
contribution to A^add becomes independent of the index 
V whereas for sufficiently small lattice spacing only j/ = 
contributes significantly. 

For small lattice spacing, the distribution p^ has a 
Gaussian shape with a width of 2a, and for a ^ 1, it 
develops a plateau with sharp edges at ±8a^, cf. Fig. [51 
The height of pr at the origin scales like 0^"+^ for small 
lattice spacing and like a~^ for large a. 

The distribution of chirality over the real eigenvalues 
p^ is shown in Fig. [5l For small a we observe the spec- 
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FIG. 4: Log-log-plot of the additional real eigenvalues versus 
a for various v. The number of matrices and its size vary 
in this plot for the Monte Carlo simulations (symbols). The 
statistical error of the numerics varies between 0.1% and 10% 
around the analytic result (solid curves). 
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In the regime of small lattice spacing, a ~ 0.1, the 
width of the distribution of the complex eigenvalues is 
given by CT = 2a^JW^/V /S whereas the spacing of the 
projection of these eigenvalues onto the imaginary axis is 
equal to A A = tt/SF. We thus have that a/ AX = 2a/ n, 
which allows us to extract a numerical value for Ws from 
lattice simulations. 

An important result is that the number of additional 
real modes is strongly suppressed for large i'. This im- 
plies that for large volumes when most configurations 
have an index ji^l > 0, additional real modes are not 
much of a problem for lattice QCD simulations with Wil- 
son fermions provided that Wga^V <C 1. 
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FIG. 5; For small lattice spacing (solid curves) the distribu- 
tion is given by the GUE (See the legend for the values 
of the index and the lattice spacing). For a 3> 1 (dashed 
curves) the shape becomes //-independent with two peaks at 
±8a^ that behave as l/\/(8a2)2 - x'^ for S > 1. 
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FIG. 6: The eigenvalue distribution pr has a Gaussian shape 
for small a (solid lines) but becomes box like with increasing 
lattice spacing (dashed curves). Notice that pr for ({v,a) — 
(0,0.1)) and ((z^, a) = (1,0.1)) is one order larger and one 
order smaller than shown in the diagram. 



tral density of the i/-dimensional Gaussian unitary en- 
semble. For large lattice spacing it deforms into a curve 
with two peaks at ±8a^ that up to an overall normaliza- 
tion is independent of ly and evolves into inverse square 
root singularities for a — !■ oo. 

Conclusions. Discretization effects become strong for 
a « 0.5. The oscillations of the spectral density in the 
continuum limit are no longer visible while the density of 
the complex eigenvalues develops a plateau with a width 
of 16a^. In terms of physical parameters, a = a^/WsV, 
with Ws a low energy constant [13] and V the volume of 
space time, we have the condition that a <C 1/^/WsV to 
be close to the continuum limit. 
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